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. ABSTRACT 
O, 

I Parton distribution functions for protons are devised in which a light gluino is 

O^l . included along with standard model quark, antiquark, and gluon constituents. 

bX)! A global analysis of a large set of hadronic scattering data provides empirical 

^ ■ constraints on the allowed range of the gluino mass as a function of the value 

I of the strong coupling strength as(Mz)- We find that m~ > 12 GeV for the 

Q>^ ■ standard model world-average value as{Mz) = 0.118. Gluino masses as small 

as 10 GeV are admissible provided that as{Mz) > 0.12, about one standard 
deviation above the world-average value. Current hadron scattering data are 
insensitive to the presence of gluinos heavier than ^ 100 — 150 GeV. 



> 

o 

1. Introduction 

(N 
00 

I Once admitted, "light" supersymmetric (SUSY) states (mass m < 100 GeV), influence 

^ ' hadron scattering processes in several ways. They change the evolution of the strong 
coupling strength as{fi) as the scale /i is varied. Second, they provide additional partonic 
p ■ degrees of freedom and share in the proton's momentum along with their standard model 
p • (SM) counterparts, altering the coupled evolution equations that govern the functional 

D ! change of the parton distributions as momentum is varied. Third, they contribute to 
^ ' hard-scattering processes as incident partons and/or as produced particles. For example, 
gluinos g materialize as hadron jets and may increase the rate for jet production at 
^ ' large values of transverse energy Et- Within the context of a global analysis of hadron 
' scattering data, these three influences allow potentially strong constraints to be placed 
on the existence and masses of SUSY particles. In this paper, I summarize the results of 
a new study of bounds on the gluino mass from a global parton density function (PDF) 
analysis p. 

In our study, we consider the effects of a light gluino only, and we do a series of 
flts to the full set of hadron scattering data used in the CTEQ6 PDF analysis j2] for 
various gluino masses rrig and values of as{Mz)- All SM contributions are evaluated in 
perturbative QCD at next-to-leading order (NLO). The g contributions to the parton 
density evolution and jet production cross section are evaluated at Born level (sufficient 
because g{x) <C q{x), g{x)). Data on hadron jet production from the Fermilab Tevatron, 
pp ^ jet -|- A at large Et are included in the fit. Quantitative bounds on m~ are obtained 
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from a PDF error analysis similar to that of CTEQ6. 

1.1. Role ofa,{Mz) 

Our conclusions on the allowed values of rrig depend on what is assumed for as(M^), 
the value at which the evolution of as{^) is pinned, but we must and do fit data at all /i. 
In our PDF analysis, we perform a series of fits to the hadron scattering data in which 
as{Mz) is varied over the broad range 0.110 < as{Mz) < 0.150. We then impose the 
direct Z-pole two standard deviation {2a) bounds on as{Mz) to establish bounds on rrtg. 

A combined analysis of all Z-pole data within the context of the SM provides as{Mz) = 
0.1187 ± 0.0027 |3j. The SM world-average obtained from a variety of determinations at 
different scales fi yields as{Mz) = 0.1183 ± 0.0027 jlj. Although we use these values to 
obtain our bounds on rrig, we are mindful that SUSY-QCD contributions can affect the 
extraction of as{fi). A recent estimate [HI finds as{Mz) = 0.118 — 0.126±0. 005, where the 
variation in the central value arises from uncertainty in the size of SUSY-QCD corrections 
to the Z width from, e.g., Z — >■ bb*g/bbg; b denotes a bottom squark. As remarked later, 
a larger value of as{Mz) permits smaller values of rrig. 

1.2. Parton distribution function (PDF) analysis 

The PDF's /(x, Q) are parametrized at a starting reference scale Qo and evolved to all 
Q > Qq. The gluino distribution is generated radiatively from gluon splitting for Q > rrig. 
The coupled evolution equations in the presence of a light gluino may be found in Ref . jT] . 
Normalized cross sections are calculated based on the parton distributions. Agreement 
with experiment is measured by x^. The PDF shape parameters aX Q = Qq are varied 
to minimize y^. The gluino parton density remains very small with respect to the SM 
quark and gluon densities for all accessible values of Q. The gluinos remove momentum 
primarily from the gluon density, for values of x > 0.05. Since hadron jet data at large 
values of Et are known to probe the gluon density at large x, inclusion of these data in 
the fit plays a strong role in constraining rrig. 

1.3. SUSY hard-scattering contributions 

Once SUSY particles are admitted as degrees of freedom, we must consider their 
impact on all hard scattering processes. For jet production in — > jet we include 
0{al) subprocesses with gluinos produced in the final state and either one or no gluinos in 
the initial state: q+q —>■ g+g, g+g —>■ g+g, g+g —>■ g+g, and q+g — > q+g. Contributions 
to the cross section from these processes are added to the SM contributions, and 

they offset partially the loss of jet rate from the smaller gluon PDF. In deep-inelastic 
scattering and in massive-lepton-pair production, the SM processes contribute at 
but the SUSY processes enter at 0{al). We can neglect 0{al) SUSY subprocesses such 



as 7* + fl' 9 + 5') where q denotes a squark, and 0{a1) subprocesses Y + q-^q + g + g- 
2. Results 

Our principal results are presented in the form of a contour plot. We show the differ- 
ence Ax^ between the value of obtained in our fit and that of a purely SM fit equivalent 
to the CTEQ6M fit. The CTEQ6 tolerance criterion for an acceptable fit is A^^ < 100; 
the isoline corresponding to Ax^ = 100 is shown in Fig.Q]by the solid line. The acceptable 
fits lie inside a valley that extends from large gluino masses and as{Mz) = 0.118 down 
to rrig ~ 0.8 GeV and right to as{Mz) = 0.145. An even narrower area corresponds to 
fits with close to those in the CTEQ6M fit. We note that is better than in the 
CTEQ6M fit in a small area in which nig < 20 GeV and as{Mz) > 0.125, with the min- 
imum Ax^ ^ —25 at nig = 8 GeV and as{Mz) = 0.130. While perhaps intriguing, this 
negative excursion in Ax^ is smaller than the tolerance Ax^ = 100 and cannot be inter- 
preted as evidence for a light gluino. We determine that m~ > 12 GeV at as{Mz) = 0.118. 
The lower limit increases as as{Mz) rises; for example, with as{Mz) = 0.122, ~ 1 a above 
the SM world average, m~ > 5 GeV. 
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Figure 1: Contour plot of Ax^ = [as{Mz),mg] — X(]'pjT;Qg]y[ obtained from our fits to hadronic 

scattering data for various nig and as{Mz)- The value of as{Mz) serves as the abscissa, with the gluino 
mass along the ordinate. The vertical dashed line marks the SM value as{Mz) — 0.118. 



3. Comparison with other bounds on m~ 

Our global analysis is based on theoretically clean one-scale observables only. We 
make no assumptions about the stability of the superpartners. Our bound m~ > 12 GeV 
is stronger than the bound nig > 6.3 GeV at as{Mz) = 0.118 found in a study of the Z 
boson width measurements 0. Negative searches for stable hadronizing gluinos provide 
bounds of nig > 18 GeV [7j and rrig > 26.9 GeV jH], but these bounds do not apply if the 
gluino decays through an i?-parity violating mechanism [Oj. Based on a study of jet shape 
variables, the DELPHI collaboration derives a limit rrig > 30—40 GeV [lOj. However, their 
analysis deals with multi-scale observables, and the treatment of theoretical uncertainties 
may be too optimistic. 

4. Conclusions 

In our work, we obtain model-independent bounds on the existence and mass of color- 
octet fermions (gluinos) based on a fit to the complete CTEQ6 set of inclusive hadron 
scattering data (1811 points). The data are characterized by high precision and cover a 
broad range in x and Q^. We include SUSY contributions to the Tevatron jet production 
cross sections and use the full CTEQ error analysis to obtain our bounds. We determine 
that rrig > 12 GeV for as{Mz) = 0.118, with smaller values of rrig allowed for larger 
as{Mz). The possibility [TT] remains open for rrig = 10 — 20 GeV if as{Mz) > 0.119. 
For as{Mz) > 0.124, rrig is also bounded from above. The PDF analysis of current 
hadron scattering data is not sensitive to gluinos with mass rrig above the weak-scale, but 
inclusion of a gluino with rrig ~ 100 GeV slightly improves the description of jet data at 
high Et pp. Our study is complementary to those in which bounds rrig are derived from 
LEP data, and it demonstrates the potential for PDF analysis to independently constrain 
new physics in the next few years once high precision data extend to larger values of Q. 
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